VERITAS Observations of the 7-ray Binary LS I +61 303 
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ABSTRACT 

LS I +61 303 is one of only a few high-mass X-ray binaries currently detected 
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at high significance in very high energy 7-rays. The system was observed over 
several orbital cycles (between September 2006 and February 2007) with the 
VERITAS array of imaging air-Cherenkov telescopes. A signal of 7-rays with 
energies above 300 GeV is found with a statistical significance of 8.4 standard 
deviations. The detected flux is measured to be strongly variable; the maximum 
flux is found during most orbital cycles at apastron. The energy spectrum for the 
period of maximum emission can be characterized by a power law with a photon 
index of T = 2.40 ± 0.16 sta t ± 0.2 sys and a flux above 300 GeV corresponding to 
15-20% of the flux from the Crab Nebula. 

Subject headings: gamma rays; LS I +61 303; X-ray binary; acceleration of 
particles 

LS I +61 303. 



1. Introduction 



The high-mass X-ray binar y LS I +61 303 consists of a massive Be-type star surrou nded 
by a dense circumstellar disk (IHutchings fc Crampton 1 1 1981c iGregory fc Neish 1 120021 1 and 
a compact object. It has rece ntly been detected as a source of very high ener gy (VHE) 7- 
rays by the MAGIC telescope ( klbert et allbood ) and confirmed by VERITAS jMaier et al. 



2007 ). This increases th e number of 7- ray binaries to three; th e other two are PSR B1259-63 
( lAharonian et al.ll2005l ) and LS 5039 (lAharonian et al.ll2006al ). While the 7-ray emission in 
PSR B1259-63 is powered by the relativistic wind of the young 48-ms pulsar, the unknown 
nature of the compact objects in LS I +61 303 and LS 5039 does not exclude microquasar- 
type emission models. These X-ray binary systems provide unique laboratories for studying 
particle acceleration by providing detailed information about the time-evolution of the par- 
ticle spectrum over the orbital period. 

Optical and radio observations of LS I +61 303 show that the compact object orbits 
the massive star every 26.496 days; the elliptical orbit is characteriz ed by a semi-major 



axis of only a few ste llar radii and an eccentricity of 0.72 ± 0.15 (ICasares et al. 



2005 



Grundstrom et al.ll2007l ). The distance to LS I +61 303 is approximately 2 kpc (jSteele et al. 
19981 ). The radio, optical and X-ray emission of the binary has a component which is 



clearly modulated at the orbital period (ITaylor &: Gregory! Il982l ; iMendelson fc Mazehl 11989 



Leahy! l200ll ; I Wen et al.l 120061 ; ISmith et al.l 12007), although consid erable variation in the 



lightcurve is observed from orbit to orbit (e.g. ISidoli et al.l (120061 )). Orbital mod ulation 
in low-energy (>100 MeV) or high-energy 7-rays (>300 GeV) is not yet confirmed (IMassi 
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20041 ; lAlbert et al.ll2006r). An associ ation of LS I +61 303 with the COS-B 7-ray s ource 
2CG 135+01 (ISwanenburg et al.lll98ll ) was proposed early on (ITaylor &: Gregoryll 19781 ). and 
obse rvations with EGRE T (3EG J0241+6103) for photon energies >100 MeV supported 
this (IKniffen et al.l 119971 ). The MAGIC observations in VHE 7-rays provided the first firm 
connection of the site of the variable very-high energy 7-ray emission with LS I +61 303. 

This paper reports on stereoscopic observations of LS I +61 303 with the ground-based 
7-ray observatory VERITAS, at energies above 300 GeV. A publication descri bing contem- 
poraneous X-ray and TeV observations is in pre paration (lAcciari et al.ll2008l ). We adopt 
orbital parameters as derived from radio data by [Gregory J2002h with an orbital period of 
P = 26.4960 ± 0.0028 days and zero orbital phase of T = HJD 2443366.775. The periastron 
takes place at phase 0.23, apastron is at phase 0.73, and infe rior and superior conjunctions 
are at phases 0.26 and 0.16 respectively (jCasares et al.ll2005l ). 



2. Observations and analysis 

The observations of LS I +61 303 were made between September 2006 and February 
2007 during the construction phase of VERITAS, an array of four 12-m imaging Cherenkov 
telescopes. VERITAS is located at the basecamp of the Fred La wrence Whippl e Obse rvatory 



in southern Arizona (1268 m a.s.L, 31°40'30"N, 110°57'07"W, IWeekes et all (12002( 1) . The 
system combines a large effective area (> 3 x 10 4 m 2 ) over a wide energy range (100 GeV to 
30 TeV) with good energy (10-20%) and angular resolution (< 0.14° on an event-by-event 
basis). 

The overall design of all four VERITAS telescopes is i dentical. Each telescope employs a 
12 m-diameter tessellated mirror of Davies-Cotton design (jDavies fc Cotton 1957 ) with 12 m 



focal length, mounted on a altitude-over-azimuth positioner. The reflector comprises 350 
hexagonal mirror facets giving a total mirror area of 106 m 2 . The focal plane is equipped 
with a 499-element photomultiplier-tube (PMT) imaging camera. The angular pixel spacing 
is 0.15°, giving a field of view of 3.5°. Light cones installed in front of the cameras increase the 
photon collection efficiency and shield the PMTs from ambient light. The three-level trigger 
system of the VERITAS array allows a substantial suppression of background events at the 
trigger level. It especially suppresses events due to local muons which are a considerable 
background in single-telescope operation. The first level of the trigger system consists of 
custom built constant fraction discriminators (CFD), one for each PMT. All observations 
described here were made with a CFD threshold of 50 mV, corresponding to approximately 
4-5 photoelectrons. The second level, a pattern trigger, requires at least three adjacent 
triggered pixels in order to generate a camera trigger. The array trigger determines if level- 
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two triggers from individual telescopes are consistent with an air shower. A coincidence 
of at least two telescopes triggering within a time window of 100 ns is required. When the 
array is triggered, PMT signals in each telescope are digitized using 500 Megasample/second 
custom-built flash-ADC electronics. To achieve a large dynamic range, an autoranging gain 
switch extends the dynamic range from 256 to 1500. In the event of a telescope trigger, all 
signals (i.e. a 48 ns long flash-ADC trace) from all channels in all telescopes are read out as 
a stereo eve nt. A description of the technical deta ils and the performance of VERITAS can 



be found in ([Holder et al.ll2006l ; iMaier et al.l 120071 ) and references therein. 



Observations of LS I +61 303 with VERITAS were made with different configurations as 
the array construction proceeded; no changes to the individual telescopes were implemented 
in this period. Observations from September 2006 to November 2006 were made with a two- 
telescope system while those for January and February 2007 were made with three telescopes. 
The typical array trigger rate was 90 events/s for the two-telescope system (with a dead time 
of 4%), and 160 events/s (7% dead time) for the three-telescope system. Data were taken on 
moonless nights in "wobble-mode" , wherein the source was positioned at a fixed offset from 
the camera center. This has the advantage that no off-source observations are necessary for 
background estimation. The conservative energy threshold after the analysis cuts adopted 
here is 300 GeV at 30° zenith angle. The dataset consists of 45.9 hours of observations 
after quality cuts and dead-time corrections. The quality cuts remove all data taken under 
non-optimal conditions (e.g. bad weather, runs with large trigger-rate fluctuations, runs 
with malfunction of the telescope system). There is good coverage of LS I +61 303 for 
orbital phases from 0.4 to 0.9, but since observations were precluded by the bright moon, 
no observations were possible during orbital phases from 0.9 to 0.2. Table [1] shows the most 
important observational parameters for each of the five observation periods. 

The first step in the analysis is the cal ibration and integration of the flash- ADC traces, 



as described in detail in (IHolder et al.ll2006l ). The calibration is divided into several sections. 
The absolute calibration uses a laser system in order to determine the signal size produced 
by single photons. In the relative calibration, laser events are used to calculate the relative 
gains of the pixels and timing differences due to path length differences in the cabling of 
each channel. Pedestal events, injected at 1 Hz during an observation run, are used for an 
estimation of the voltage offset in each flash-ADC trace and the noise levels due to the night 
sky background and electronics. The amount of charge in each flash-ADC trace is calculated 
by summing the samples for a given window size. A two-pass summation method is used 
here, allowing the integration of the flash-ADC traces with an optimal signal-to-noise ratio. 
The resulting image of an air shower is cleaned in order to remove pixels which contain 
mainly background light. The cleaning consists of a two- level filter removing all pixels with 
an integrated charge smaller than 5 times their pedestal standard deviation and any pixel 
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that are adjacent to these higher threshold pixels and having signals smaller than 2.5 times 
their pedestal sta ndard devia tion. The shower image is then parameterized with a second 
moment analysis (IHillas Ill985l ). The direction of origin of the 7-ray on the sky and the impact 
parameter of the shower core on the ground are reconstructed using stereoscopic techniques 
(IHofmann et al.l Il999l ; iKrawczynski et al.l 120061 ) . At least two images with an integrated 
charge per image > 400 digital counts (~75 p.e.) and a maximum image distance from the 
center of the camera of less than 1.2° are required in this reconstruction stage. The majority 
of the far more numerous background events are rejected by comparing the shape (i.e. width 
and length) of the event images in each telescope with the expected shapes of 7-ray showers 
modeled by Monte Carlo simulations . These so-called mean-s caled width and mean-scaled 
length parameters fjKonopelko I Il995l ; IKrawczynski et al.l 120061 ) are calculated with lookup 
tables based on Monte Carlo simulations. The lookup tables contain the median and 90%- 
widths of the image parameter width (w M c, Vwidth,Mc) and length (l MC , cri en gth,Mc) as a 
function of impact parameter R, integrated charge per image s, and zenith angle G: 
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and similar for mean scaled length. The cuts applied here are —1.2 < mean scaled width/length 
< 0.5 and reconstructed distance of shower core position from the center of the array < 
250 m. This and an additional cut on the arrival direction of the incoming 7-ray (@ 2 <0.025 deg 2 ) 
reject more than 99.9% of the cosmic-ray background while keeping 45% of the 7-rays. All 
the cuts are optimized a priori with Monte Carlo simulations of 7-ray- and hadron-induced air 
showers. The background in the source region is estimated from the same field o f view using 



the "r eflected-region" model and the "ring-background" model as described in iBerge et al. 
j2007h . 



The energy of each event in the source and background region is estimated from Monte 
Carlo simulations assuming that the primary particle is a 7-ray. The calculation uses lookup 
tables and determines the energy of an event as a function of impact parameter, integrated 
charge per image, and zenith angle. Collection areas at different zen i th an gles for 7-rays are 
produced from Monte Carlo simulations following (IMohanty et al.lll998l ); collection areas 
are interpolated between zenith angle bins to the zenith angle of each event. The limited 
energy resolution is taken into account by calculating the collection areas as a function of 
reconstructed energy. The differential flux dN/dE in each energy bin of width AE is then 
calculated with the number of events in the source region (N source ) and in the background 
region (Nb C k), the ratio of the size of the source region to the size of the background region 
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(a), the dead time corrected observation time T, and the collection area A: 

(Nsource N bck \ 
^ 1/A - or 1/A /T/AE 
i=i i=i / 

The dependence of the collection area on the spectral index is taken into account by an 
iterative process. In each step of the iteration, collection areas are calculated using the spec- 
tral index obtained in the previous step. The iteration stops when convergence is achieved, 
usually after 2-3 steps. In order to accurately calculate source fluxes and energy spectra, a 
model of the telescope response to air showers has been developed. The Monte Carlo simula- 
tions take all relevant processes and efficiencies in the development of the air shower through 
the atmosphere, the propagation of Cherenkov photons through the optical s ystem of the 



telesc opes, and the response of the camera and electronics into account; see ({Holder et al. 



20061 ) for more details. The systematic error in the energy estimation is dominated by un- 
certainties and variabilities of the atmospheric conditions, in the Monte Carlo simulations, 
and in the overall photon collection efficiency. 



3. Results 

We have detected LS I +61 303 as a source of 7-rays with energies above 300 GeV, 
at a total significance of 8.4 a at the position of the optical counterpart. The source has 
been found to be variable, as earlier measurements by MAGIC suggest. The probability 
that the measured fluxes from LS I +61 303 are constant with time or orbital phase has 
been determined with a x 2_ t es t to be less than ~ 10 -9 . The two-dimensional sky map of 
significances for the region around LS I +61 303 shows a strong detection for orbital phases 
0.5-0.8, i.e., around apastron (Figured]). No signal has been found for orbital phases 0.8-0.5. 
The position of the peak of the 7-ray emission, reconstructed by a fit to the uncorrelated map 
of excess events with a two-dimensional normal distribution, is in agreement with the position 
of the optical source: A R a = 147±73 sta t arcsec, A dec = — 34±30 sta t arcsec (corresponding to 
RAj 20 oo 2 h 40 m 41 s , Dec J20 oo 61°13'12"). The systematic uncertainty is ±90 arcsecondfl This 
is consistent with the position of the TeV excess reported by MAGIC. The morphology of the 
excess is compatible with the expected distribution due to a point source. The VERTIAS 
source name is VER J0240+612. 

In Figure [2]we show the differential energy spectrum during the high-flux phases 0.5-0.8 



1 There were no optical pointing monitors installed during these observations with the array under 
construction. The pointing accuracy is expected to improve to <15 arcsec for the completed system. 
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for 7-ray energies between 300 GeV and 5 TeV. The integrated flux of VHE photons in 
this energy range is (8.13 ± 1.02) x 10~ 12 cm~ 2 s _1 . The shape is consistent with a power 
law dN/dE = Cx(E/l TeV)~ r with a photon index T = 2 AO ± 0.19 stat ± 0.2 sys and a flux 
normalization constant C = (2.89 ± 0.32 stat ± O.Q sys ) x 10~ 12 cm-^TeV- 1 . The x 2 of the 
fit is 2.1 for 5 degrees of freedom. 

The dependency of integral fluxes on orbital phase (Figure [3]) shows that significant 
fluxes above 2a are only detected during orbital phases close to apastron. Peak fluxes 
measured at phases between 0.6 and 0.8 correspond to 10-20% of the flux of the Crab 
Nebula, but appear to vary from one orbital cycle to the next. Upper limits in a range 
corresponding to 3-10% of the flux of the Crab Nebula have been derived for phases 0.3 
to 0.6 and around 0.9. The upper limit calculation assumes a similar shape for the energy 
spectra of TeV 7-rays during all phases. This may not be correct , as the phase- dependen t 



PC 

2 



shape of the energy spectra of the X-ray binary LS 5039 shows (lAharonian et al.l l2006bl ). 
Furthermore, the combination of data from different orbital cycles in Figure [3] (bottom) 
assumes implicitly that the position of the emission maximum in orbital phase and its flux 
do not vary from one cycle to the next. A constant maximum 7-ray flux is not necessarily 
expected; LS I +61 303 shows variable X-ray fluxes over different orbits (e.g. Chernyakova 
et al. 2006). The VERITAS measurement can neither confirm nor refute this at this stage. 
The same periodicity for radio and 7-ray data was assumed during this a nalysis (P=26.4960 



days). This has been tested using Lomb-Scargle statistics (IScargldll982l ). but large gaps in 
data taking due to bright moon periods and the periodicity of the observation cycles do not 
allow definite conclusions. 



Discussion and Conclusions 



LS I +61 303 has been detected in 7-rays at phases around 0.6 to 0.8 only, when the 
distance between the two objects in the binary system is largest. This indicates a strong 
dependence of particle acceleration and energy loss mechanisms on the relative positions. 

It is suggested that TeV 7-rays are produced in LS I +61 303 in leptonic (inverse 
Co mpton scattering of low-energy photons fro m the stellar companion or pair ca scades, see 



Bednarekl (120061 ) ; lGupta fc Bottcherl (120061 )) and/or hadronic interactions (e.g. lRomero et al 



( 120051 )). Both production mechanisms require a population of high-energy particles and a suf- 
ficiently dense photon field or target material. This is provided abundantly in LS I +61 303 
for all orbital phases by stellar photons and the wind of the companion star. Target density 
is highest around periastron, which favors a maximum in 7-ray emission when the compact 
object passes close to the stellar companion. The maximum electron energy available for 
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7-ray production is, on the other hand, restricted by cooling through inverse Compton and 
synchrotron emission; both of these effects can be strongest around periastron. In addi- 
tion, 7-rays emitted too close to the massive star suffer from photon-photon absorption in 
the dense stellar radiation field. Significant i nhomogeneities in the stellar w i nd can further 



compl icate the process of 7-ray production (INeronov fc Chernyakoval 120071 ; iRomero et al. 
2007bh . 



The unknown nature of the compact object (neutron star or black hole) allows two very 
different scenarios to provide the necessary power for the production of VHE 7-rays. In the 
micro quasar model, c harged particles (electrons or hadrons) are accelerat ed in an accretion- 
driven relativistic jet (ITaylor &: Gregorylll984l ; iMirabel &: Rodriguezlll994t). Acceleration can 



take place all along the jet or in shocks created in the jet termination zone (IHeinz fc Sunyaev 



20021 ). The efficiency of particle acceleration and 7-ray production varies over the orbit 
due to changes in target density, accretion rate, and magnetic field strength. The second 
scenario assumes that the compact object is a pulsar. Particles are accelerated in the shock 
created by the c ollision of the expanding p ulsar wind with the equatorial disk or wind of the 
companion star (iMaraschi fc Treves! 1 198 ll ). similar to the binary system PSR B1259-63. Th e 
distance of the pulsar wind termination shock changes with orbital motion (IDubusl 120061 ) . 
Electron cooling through synchrotron emission is in this model lowest around apastron, 
when the shock is furthest away from the pulsar and therefore the magnetic field is weakest. 
This results, under the assumption that enough target photons are available for inverse 
Compton scattering, in the observed maximum of the 7-ray emission at apastron. More 
measurements, especially contemporaneous multiwavelength observations are necessary for 
constraining tests of the available models. 

The measured integral flux of 7-rays above 300 GeV in the phase interval 0.6 to 0.8 
corresponds to an isotropic luminosity of roughly 10 34 ergs/s at 2 kpc. The necessary power 
could either be provided by a pulsar with spindown luminosity of ~ 10 36 erg/s or by 
accret ion, which has been estimated to be in the range of several 10 37 erg/s (IRomero et al. 
2007al ). Both scenarios imply conversion factors to 7-ray emission below 1%. 



If we compare the average TeV flux with the average GeV flux measured by EGRET 
( IKniffen et al.lll997l ). we find a GeV/TeV photon flux ratio of 10 5 which is consistent with 
that produced by a power-law with photon index 2.4, similar to the photon indexes measured 
in the two bands. This result should be treated with caution since the source is known to 



exhibit long-term variability ([Gregory! Il989l ) and the measurements are separated by more 
than 10 years. However, it may suggest that a single power-law could suffice to describe the 
GeV/TeV spectrum - a possibility that can be tested with future simultaneous GeV/TeV 
observations. 
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In conclusion, new stereoscopic observations with VERITAS have confirmed very-high- 
energy 7-ray emission from LS I +61 303. The source is detected only during orbital phases 
0.6 to 0.8, which is close to apastron. This suggests variability connected to the orbital move- 
ment, but a definite conclusion has to wait for further measurements. Future observations 
with the much more sensitive full VERITAS array will allow us to study the VHE emission 
during the whole orbital cycle. 
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tional Science Foundation, the Smithsonian Institution, by NSERC in Canada, by Science 
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Table I. Details of the VERITAS observations of LS I +61 303. 



Month 


N T ei a 


Wobble 


Elevation 


Obs. 


Orbital 






offset 


range 


time b 


phase 










[h] 


interval 


09/2006 


2 


0.3° 


53-61° 


12.5 


0.31-0.69 


10/2006 


2 


0.3° 


61-61° 


7.2 


0.51-0.78 


11/2006 


2 


0.3° 


54-61° 


12.8 


0.23-0.87 


01/2007 


3 


0.5° 


55-61° 


9.9 


0.41-0.76 


02/2007 


3 


0.5° 


52-57° 


3.5 


0.51-0.89 


Total 






52-61° 


45.9 


0.23-0.89 



a Number of available telescopes 
b Dead-time-corrected observation time 
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Fig. 1. — Significance map of the region around LS I +61 303 in equatorial J2000 coordi- 
nates. Left: Observations during orbital phases 0.8 to 0.5 (18.5 h of data). Right: Observa- 
tions during orbital phases 0.5 to 0.8 (around apastron, 25 h of da ta). Significances are n ot 
corrected for number of trials. The position of the optical source (jPerryman et al. 1997 ) is 
indicated by a white cross. The background is estimated using the ring-backgrou nd model- 
Neigh boring: bins are correlated. Significances are calculated using the method of iLi fc Ma 
(119831 ). equation 17. 
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Fig. 2. — Differential energy spectrum of VHE photons above 300 GeV for LS I +61 303 
around apastron (orbital phases 0.5-0.8). The markers indicate measured data points, the 
continuous line a fit assuming a power-law distri bution of the data. Downward pointing 
arrows indicate upper flux limits (99% probability. iHelend (119831 )) for bins with significances 
below 2a. For comparison, the energy d istributions of LS I +61 303 published by the MAGIC 
collaboration for orbital phases 0.4-0.7 ( Albert et al. 20061 ) and of the Crab Nebula measured 
by VERITAS (September to November 2006) (reconstructed as dN/dE = 2.85 x 10~ u ■ 
E _2 - 49 cm _2 s~ 1 TeV _1 ) are indicated by dotted and dashed lines, respectively. 
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Fig. 3. — Average fluxes or upper flux limits per orbital phase bin for 7-rays with energies 
above 300 GeV from the direction of LS I +61 303 as a function of orbital phase. The 
bottom panel shows the results averaged over the whole dataset, the upper pan els show the 
results for individual orbits. Upper flux limits (95% probability, Helene ( 1983h N ) are show n 
for data points with significances less than 2a (significance calculation after iLi fe Mai (119831 ). 
equation 17). A flux corresponding to 10% of the flux of 7-rays from the Crab Nebula is 
indicated by the dashed horizontal line in the bottom panel. 



